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ABSTRACT 

We announce the discovery of a new eclipsing hot subdwarf B + M dwarf binary, 
EC 10246-2707, and present multi-colour photometric and spectroscopic observations 
of this system. Similar to other HW Vir-type binaries, the light curve shows both 
primary and secondary eclipses, along with a strong reflection effect from the M dwarf; 
no intrinsic light contribution is detected from the cool companion. The orbital period 
is 0.118 507993 6 ± 0.000000 0009 days, or about three hours. Analysis of our time- 
series spectroscopy reveals a velocity semi- amplitude of Ki = 71.6 ± 1.7 km s _1 for 
the sdB and best-fitting atmospheric parameters of T c ff = 28900 ± 500 K, log g = 
5.64 ± 0.06, and log N(He)/N(H) = -2.5 ± 0.2. Although we cannot claim a unique 
solution from modeling the light curve, the best-fitting model has an sdB mass of 0.45 
Mq and a cool companion mass of 0.12 Mq. These results are roughly consistent with 
a canonical-mass sdB and M dwarf separated by a ~ 0.84 Rq. We find no evidence of 
pulsations in the light curve and limit the amplitude of rapid photometric oscillations 
to < 0.08%. Using 15 years of eclipse timings, we construct an O-C diagram but find 
no statistically significant period changes; we rule out \P\ > 7.2 x 10~ 12 . If EC 10246- 
2707 evolves into a cataclysmic variable, its period should fall below the famous CV 
period gap. 
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1 INTRODUCTION 

Few stars seem to prefer companionship as much as the hot 
subdwarf B (sdB) stars. It is widely accepted that these 
evolved, post-main sequence stars are the progeny of red 
giants that were somehow stripped of their outer hydrogen 
envelopes, leaving behind a 0.5 Mq helium-burning core sur- 
rounded by a thin layer of hydrogen (Heber 1986). Mengel 
et al. (1976) first proposed binary interactions as a probable 
cause for the stripping. Since then, various models have been 
constructed around this idea, most of which use the angular 
momentum stored in orbits to spin up and eject the progen- 
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itor's envelope. Observations tend to support this picture, 
with reports of the sdB binary fraction ranging from 20% to 
nearly 100% over the past few decades (see Heber 2009 for 
a review). Today, the exact value remains uncertain, but at 
least a small fraction of sdBs are thought to be single stars. 

The Han et al. (2003) models highlight five main for- 
mation scenarios leading to subdwarf B stars: the "first" 
and "second" Roche lobe overflow channels, the "first" and 
"second" common envelope channels, and the merger of two 
helium white dwarfs. More recently, Clausen & Wade (2011) 
proposed a channel in which a helium white dwarf and a low- 
mass, hydrogen-burning star merge to create a hot subdwarf. 
Since binary population synthesis models predict specific 
distributions of orbital periods, companion types, and sub- 
dwarf masses, they may be tested by observations and mea- 
surements of these parameters. Unfortunately, the masses of 
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sdBs are difficult to derive; fewer than 15 have been deter- 
mined to date. Most were derived from asteroseismological 
analyses, and, except in the case of PG 1336-018 and 2M 
1938+1946, they have not been tested using independent 
techniques. Other techniques leading to masses include the 
analysis of double-lined spectroscopic binaries, light-curve 
modeling of eclipsing binaries, and transit timings of eclips- 
ing binaries. 

This work focuses on the short-period sdB binaries, 
which are thought to be products of common envelope evo- 
lution. Dynamically unstable mass transfer between the hot 
subdwarf progenitor, a red giant, and its lower-mass main 
sequence companion leads to the formation of a common en- 
velope and, eventually, the expulsion of the red giant's outer 
hydrogen envelope. This process produces a hot subdwarf 
with a cool main sequence companion in a tight, circular or- 
bit. Orbital parameters have been measured for more than 
100 close binaries and show periods ranging from 0.07 to 28 
days with a median near 0.61 days (see Table A.l of Geier 
et al. 2011; Copperwheat et al. 2011). Most measured peri- 
ods fall below 1 day; only a handful are longer than 10 days. 
We refer the reader to Figure 5 of Barlow et al. (2012) for 
an up-to-date orbital period histogram of all sdB binaries 
with solved orbits. 

Around a dozen sdB binaries are currently known to 
show eclipses; most have M dwarf companions. Although 
the cool companions contribute a negligible amount to the 
total system flux, these binaries show both primary eclipses 
(when the M dwarf occults the sdB) and secondary eclipses 
(when the sdB occults the M dwarf) ; in the latter case it is 
the reflection effect from the companion which is occulted. 
It is widely believed these systems will evolve into catac- 
slysmic variables (CVs) once the orbits decay sufficiently 
from angular momentum loss (see Warner 2003 for a re- 
view). Qian et al. (2008), however, present evidence show- 
ing the orbits might evolve so quickly that mass transfer 
begins before the hot subdwarf has evolved into a white 
dwarf, leading to something other than a cataclysmic vari- 
able. Long-term monitoring of eclipse timings will help shed 
light on their evolutionary rates. Since modeling the binary 
light curves can lead to estimates of the component masses, 
each additional eclipsing sdB system discovered helps im- 
prove our understanding of the sdB mass distribution of 
common envelope-produced binaries. 

Here, we present estimates of the system parameters 
for the new eclipsing sdB+dM system EC 10246-2707 (here- 
after, EC 10246; a = 10 ,1 26 m 56.59 s , S = -27°22'58.73", 
J2000; B=14.2), which was discovered to be an sdB star in 
the Edinburgh-Cape survey (Kilkenny et al. 1997). Two of 
us (BNB &: BHD) discovered photometric variations in EC 
10246 on 21 May 2009 while carrying out a survey for new 
pulsating sdB stars with the 0.4-m Panchromatic Robotic 
Optical Monitoring and Polarimetry Telescopes (PROMPT) 
on Cerro Tololo in Chile (Reichart et al. 2005). Unbeknownst 
to them at the time, two of the other co-authors had al- 
ready discovered eclipses at SAAO (DO'D) and had been 
observing the system over several seasons (DO'D and DK). 
Using simultaneous spectroscopy and multi-colour photome- 
try, we derive light-curve modeling solutions for the system 
and draw comparisons to other known eclipsing sdB+dM 
binaries. 

In §2, we present time-series spectroscopic observations 
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UT Date 


15 Dec 2009 


Start Time a 
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Run Length 


3.1 hr 
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Slit Width 
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Dispersion 


0.65 A pixel- 1 
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40 s 




Dead Time 


9 s 




Cycle Time 


49 s 




Duty Cycle 


82% 




No. exposures 


228 





"mid-exposure time 



and the radial velocities and atmospheric parameters de- 
rived from them. Photometry of the system is described 
in §3, along with an orbital ephemeris and the details of 
our binary light curve modeling. §4 follows with a presenta- 
tion of the best-fitting system parameters derived from this 
modeling. Finally, we discuss several miscellaneous items in 
§5, including similarities with other eclipsing sdB+dM bi- 
naries, orbital period changes, the subsequent evolution of 
EC 10246, and a newer method for obtaining independently- 
derived masses of the individual components. We conclude 
with a brief summary of our work in §6. 

2 TIME-RESOLVED SPECTROSCOPY 
2.1 Observations and reductions 

We obtained 228 low-resolution spectra of EC 10246 with 
the Goodman Spectrograph (Clemens et al. 2004) on the 
4.1-m SOuthern Astrophysical Research (SOAR) telescope 
on 15 December 2009. We used a 10" slit with the hopes of 
creating a light curve from our spectroscopy. With such a 
large slit, the spectral resolution is set not by the slit width 
but by the seeing, which ranged from 1.0-1.5". These con- 
ditions in combination with the 600 mm -1 VPH grating 
(1.3 A per binned pixel dispersion) gave us an average spec- 
tral resolution of 5.4 A over the wavelength range 3600-6276 
A. The slit axis was aligned to 216.6 deg E of N, allowing 
us to place a second star on the slit (1.22' away) to mon- 
itor drifts in the wavelength solution over the run length. 
We read out only a subsection of the CCD (approximately 
2071x623 binned pixels) to minimize the time between ex- 
posures. Each spectrum was integrated for 40 s, resulting in 
a signal-to- noise (S/N) ratio of approximately 20 per reso- 
lution element. The entire series lasted 3.1 hours, covering 
slightly more than one orbital cycle. Table 1 presents a sum- 
mary of the observational setup. HgAr comparison spectra 
were obtained before and after the series. Using the same 
instrumental setup, we also took four 30-s spectra of the 
spectrophotometric standard star EG 21 for flux calibration. 

We used the cedproc routine in iraf to bias-subtract 
and flat-field all spectral images and apall to optimally ex- 
tract one-dimensional spectra from the frames and subtract 
a fit to the sky background. The HgAr frames we obtained 
were insufficient for wavelength calibration since the emis- 
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Figure 1. Top panel: Radial velocity curve for EC 10246, plotted 
twice for better visualization. Our full velocity curve covers 1.1 
orbital cycles, but to simplify the presentation of the data, we 
only show one continuous cycle here. Velocities were measured 
from the hydrogen Balmer profiles in the sdB spectra. The solid 
line denotes the best-fit sine wave to the data. Bottom panel: 
Residuals after subtracting the best-fit sine wave from the data. 
The mean noise level in the Fourier transform of the residual 
velocity curve is 2.2 km s _1 . 



sion lines were too broad (from the big slit) to accurately 
centroid. Instead, we used a dispersion solution from a dif- 
ferent night created from a HgAr lamp observed through 
a smaller slit. Unfortunately, this makes it impossible to 
measure absolute space velocities, although relative mea- 
surements can still be made. We applied the same reduction 
techniques to the EG 21 spectra and used them to flux- 
calibrate the EC 10246 data. The spectrum is dominated by 
hydrogen Balmer absorption features and also shows several 
helium I lines (4026 A, 4471 A, 4922 A, 5876 A). An sdOB 
classification is ruled out due to the absence of the helium 
II line at 4686 A. 



2.2 Radial velocity curve 

Relative velocity shifts were computed from the cores of 
the hydrogen Balmer profiles. We used the MPFIT routine 
in IDL (Markwardt 2009), which employs the Levenberg- 
Marquardt method, to fit inverse Gaussians to H/3-H9. The 
helium I lines were too weak in the individual spectra to 
be used for this purpose. Since the telescope guide star 
was significantly redder than EC 10246, our slit alignment 
was not constant over the observing run, resulting in a 
gradual, colour-dependent velocity shift. By approximating 
the colour of the guide star, we were able to remove this 
wavelength-dependent shift. Velocity drifts due to instru- 
mental flexure were removed by tracking the absorption-line 
features in the companion star on the slit. Figure 1 shows 
the resulting RV curve for the sdB, plotted twice for better 
visualization. The larger scatter near the peak of the curve 
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T = 28900 K 
logg =5.64 
log N(He)/N(H) = -2.5 
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Figure 2. Normalised spectrum of EC 10246 (black dot-dashed 
line) with the best-fitting model (solid red line). Absorption line 
profiles arc offset by 0.25 from one another. 



resulted from diminishing weather conditions towards the 
end of our observing run. 

We determined the semi- amplitude of the velocity vari- 
ation (K\) by fitting a sine wave to the data with MPFIT 
and fixing the orbital period and phase to the values deter- 
mined from the eclipse timings (see §3.2). We derive a line- 
of-sight radial velocity semi-amplitude of K\ = 71.6 ± 1.7 
km s _1 for the sdB component. Residuals from subtracting 
this fit, shown in the bottom panel of Figure 1, are con- 
sistent with random noise. We also fitted eccentric models 
to the data (with low eccentricity), but currently have no 
reason to prefer them over the circular model. 



2.3 Atmospheric parameters 

We fitted grids of metal-line blanketed LTE model atmo- 
spheres with solar metallicity (Heber et al. 2000) to the spec- 
tra with the Spectrum Plotting and Analysis Suite (SPAS; 
Hirsch 2009) to determine the effective temperature, surface 
gravity, and helium abundance of the hot subdwarf compo- 
nent. As some sdB+dM binaries have spectroscopic param- 
eters that reportedly vary over each orbital cycle (see for 
example, Heber et al. 2004; For et al. 2010), we analysed 
the data according to their orbital phase. We divided the 
time series into 10 phase bins, averaged the spectra in each 
bin, and fitted atmospheric models to the resulting spec- 
tra. Hydrogen Balmer lines H/3-H10 and the 4471 A helium 
I line were fitted with a \ 2 minimization technique. The 
best-fitting atmospheric parameters showed no change with 
orbital phase. Overall, we adopt the weighted average of 
their values and report T cff = 28900 ± 500 K, log g = 5.64 
± 0.06, and log AT(He)/JV(H) = -2.5 ± 0.2. In Figure 2, 
we show the best-fitting solution with the average observed 
spectrum over one of the 10 phase bins. The best-fitting sdB 
temperature cited above is used as a fixed parameter in the 
light-curve modeling discussed in §3.3. 
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3 PHOTOMETRY 

3.1 Observations and reductions 

Simultaneously with the SOAR spectroscopy, we used 
PROMPT 1, 3, 4, & 5 to obtain approximately 3.7 hrs of 
time-series photometry with the Johnson BVRI filters. Al- 
though these telescopes are primarily dedicated to GRBs 
and science education, they are eminently useful for studies 
of rapidly pulsating stars and eclipsing binary systems such 
as EC 10246. Cycle times ranged from 40 s (B) to 20 s (I) 
and duty cycles from 90% (I) to 80% (B). On the follow- 
ing night, we took three hours of photometry through the 
B filter with SOAR/Goodman using 15-s exposures and an 
81% duty cycle. Together, these five datasets represent the 
photometry used in our light-curve modeling. We obtained 
numerous other light curves with both PROMPT and the 
SAAO 1.0-m telescope for the purpose of constructing the 
ephemeris and O-C diagram shown in §3.2. 

Reduction of the SAAO CCD frames was performed 
on-line, enabling the observer to judge the quality of the 
observations and to select suitable stars as local compar- 
isons (to correct for small transparency variations). Con- 
ventional procedures (bias subtraction, flat field correction, 
etc.) were followed, with magnitude extraction being based 
on the DoPHOT program described by Schechter, Mateo, & 
Saha (1993). 

All PROMPT frames were bias-subtracted and flat- 
fielded using the ccdproc procedure in iraf. Due to non- 
negligible dark current in these data, we also dark- 
subtracted the frames in iraf. We extracted photometry 
with the hsp^nd 1 routine by choosing apertures that maxi- 
mized the signal-to-noise ratio in the light curves. Sky annuli 
were drawn around the stellar apertures to keep track of and 
remove sky counts. We also extracted photometry of nearby, 
bright comparison stars to remove sky transparency varia- 
tions and normalised the target light curves with parabo- 
las to take out atmospheric extinction effects. The resulting 
differential light curves from SOAR and PROMPT used to 
model the binary are shown in Figures 5 and 6, respectively. 

3.2 Orbital ephemeris 

In order to compute an orbital ephemeris for EC 10246, 
times of primary minima were determined using two differ- 
ent techniques. For the PROMPT eclipses, we fitted inverse 
Gaussians to the eclipse profiles using the MPFIT routine in 
IDL. The SAAO minima were determined by hand via the 
bisected chords technique (e.g., Kilkenny et al. 2000), which 
measures the mid-points of chords joining ingress and egress 
curves. Table 2 presents each measured time of minimum 
as a Barycentric Julian Date of the Barycentric Dynamical 
Time (BJDtdb; Eastman et al. 2010). From these measure- 
ments, which span more than a decade, we constructed a 
linear ephemeris of the form 

C = To + P x E (1) 

where To is a reference eclipse time, Pq the orbital period 
(at To), and E the orbital cycle number, as measured from 

1 written by Antonio Kanaan 



Table 2. Primary Eclipse Times of Minima 



Time of Minimum 


Telescope 


O-C Notes 


(BJD TDB ) 




(<0 



a combined point from PROMPT multi-colour photometry 



To (when E=0). We assessed the quality of this model by 
calculating an 'observed minus calculated' (O-C) diagram 
(see Sterken 2005; Kepler 1993). Such a plot allows us to 
place limits on the orbital period change (P) and to look for 
the presence of additional bodies in the system. 

Figure 3 shows the full O-C diagram, from which we 
derive an orbital ephemeris defined by 
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Figure 3. O-C diagram constructed from primary eclipse timings 
from SAAO (black circles) and PROMPT (blue triangles). The 
solid and dotted lines denote the best-fitting linear and quadratic 
curves, respectively. There is currently no reason to prefer a model 
with non-zero P over a constant-period model. We limit secular 
period changes to \P\ < 10 — 12 s s — 1 . 



BJDtdb = (2 455 680.562 160 ± 0.000 016) 

+ (0.118 507 993 6 ±0.000 000 000 9) x E 

We fitted models with both constant and changing periods, 
but computation of the F-statistic shows we have no reason 
to prefer the best-fitting P model over a constant -period 
solution. Accordingly, we claim neither a measurement of P 
nor the detection of tertiary companions from reflex motion. 
We can, however, limit secular period changes in the system 
to \P\ < 7.2 x 10 -12 s s^ 1 . Assuming a reasonable distance 
for EC 10246 (from the parameters derived in §3.3) and com- 
bining this with proper motion measurements from Roeser 
et al. (2010), we expect a period increase due to proper mo- 
tion on the order of 10~ 15 s s _1 (Shklovskii 1970; Pajdosz 
1995), well below our detection threshold. 

3.3 Binary light curve modeling 

Since only spectral features from the primary component 
have currently been identified, we rely on a light curve anal- 
ysis to compute the mass ratio (q) of the binary. To pre- 
pare the data for modeling, we normalised the light curves 
by the mean flux and used the calculated ephemeris to 
convert times to orbital phase. We use the standard con- 
vention of defining the primary eclipse centre as the zero 
point in phase. For the light-curve analysis, we employed 
the Modified ROche (MORO) code, which is based on the 
Wilson- Devinney code (Wilson & Devinney 1971). We refer 
the reader to Drechsel et al. (1995) for additional details 
concerning MORO. We computed separate solutions for the 
SOAR B curve ("Solution 1"), the PROMPT B curve ("So- 
lution 2"), and the simultaneous PROMPT BVRI photom- 
etry ("Solution 3"). In all cases, we fixed the binary type to 
'detached' in the Wilson-Devinney code ("mode 2"). 

We use several assumptions, boundary conditions, and 
input from spectroscopy to constrain the parameter space 
searched by the models. First, we assume the orbit is circu- 
lar, a seemingly reasonable assumption given the common- 
envelope histories of sdB+dM binaries and their purported 
short circularization timescales (Tassoul & Tassoul 1992). 
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Figure 4. Mass-radius diagrams for the secondary {top panel) 
and primary {bottom panel) stars, which illustrate the degener- 
acy inherent to the MORO light curve solutions. Each pair of 
points represents one of the 27 light— curve modeling solutions 
computed from MORO. They are colored according to solution 
number to aid pair identification in the two plots; their sizes cor- 
respond to the x 2 of each model fit to the data. The small red 
point denotes the best-fitting model. In the cool companion plot, 
the dot-dashed line represent the theoretical mass-radius relation 
for single M dwarfs from Baraffe et al. (1998), while the dashed 
line represents the mass-radius relationship derived by Knigge et 
al. (2011) for lower-main sequence stars in CVs. In the primary 
(sdB) plot, dashed and dot-dashed lines represent the 'canonical' 
sdB mass (~ 0.48 Mq) and minimum mass for He-burning, re- 
spectively. The shaded region marks our measured value of log g 
with 1-<t error bars. 



We also assume the rotation of the stars is synchronised with 
the orbit, a process models show takes only decades 2 (Zahn 
1977). Since a light-curve analysis only provides fractional 
component luminosities, we fixed the primary temperature 
(Ti) to the spectroscopic value determined in §2.3 and ad- 
justed only the secondary temperature (Th). As the radiation 
pressure in the M dwarf is insignificant, we set the secondary 
radiation pressure parameter (62) to zero. The primary lin- 
ear limb darkening coefficients (xi) for the BVRI curve were 
set to the values derived by Wade & Rucinski (1985). No 
significant evidence of a third light contribution {I3) was 
found during initial investigations of the model fitting, so 
we kept this value locked to zero when computing final so- 
lutions. Gravity darkening exponents of 1.0 and 0.32 were 
assumed for the primary (von Zeipel 1924) and secondary 
(Lucy 1967), respectively. Finally, the large reflection effect 



2 Recent observations of some sdB+dM systems (Pablo et al. 
2012, for example) imply the synchronization timescales might 
actually be longer than the EHB lifetime. 
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Figure 5. Top panel: SOAR B light curve with the best-fitting theoretical light curve, as determined from the MORO code. Only one 
full orbital cycle of the light curve is shown. Bottom panel: Residuals after subtracting the model from the data. 
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Figure 6. Simultaneous, multi-colour photometry of EC 10246- 
2707 obtained with PROMPT on 2009-12-15 (top panel). Light 
curves were taken through the B,V,R,, and I filters with PROMPT 
3, 1, 4, and 5, respectively. Although each data set covered ap- 
proximately 1.3 orbital cycles, we only show one complete cycle 
here. The best-fitting binary light curves from MORO (Solution 
#3) are plotted with solid black lines. Residuals after subtraction 
of the models are also shown (bottom panel). 



shows the secondary is cool while the lack of any other op- 
tical light from it implies it is small; these traits are con- 
sistent with a low-mass M dwarf companion. Accordingly, 
we used the theoretical mass-radius relationship for lower- 
main sequence stars from Baraffe et al. (1998) to further 
constrain our model outputs. We accept only those MORO 
solutions with mass-radius combinations for the cool com- 



panion falling within conceivable boundaries obtained from 
Baraffe et al. (1998). 

Other adjustable parameters for the modeling included 
the inclination angle (i), secondary albedo (A2), Roche sur- 
face potentials (fii, H2), mass ratio (q), primary radiation 
parameter (<5i), colour-dependent luminosities (Li, L2), and 
the secondary linear limb darkening coefficient (£1). 

MORO encountered no serious problems in the fitting 
except those connected with the usual parameter correla- 
tions. Different combinations of the inclination angle, limb- 
darkening coefficients, and ratio of the stellar radii, can pro- 
duce essentially indistinguishable light curves of nearly the 
same shape and fit quality. Thus, the code finds reason- 
able light curve solutions for a broad range of mass ratios 
(g=M2/Mi) and corresponding gravitational equipotentials 
(fli, O2). Figure 4 illustrates this degeneracy by present- 
ing primary and secondary mass-radius diagrams for the 27 
separate solutions MORO found for the SOAR B light curve 
("Solution #1"). 

The x -scm ared values of the models in Figure 4 are 
comparable, and consequently, we have no strong statistical 
reason to choose one model over the others without applying 
additional constraints. For reference, we overplot the empir- 
ical M-R relationship for lower-main sequence stars in CVs 3 
from Knigge et al. (2011), the theoretical M-R relation for 
single M-dwarfs from Baraffe et al. (1998), the 'canonical' 
sdB mass, and the minimum mass needed for He-burning. 
In the sdB panel, we also highlight combinations of masses 
and radii consistent with our spectroscopic measurement of 
the surface gravity. Several of the models fall within this 
shaded region, so the addition of our log g estimate alone is 
not sufficient for choosing a single solution. In light of the 
degeneracy inherent to the models, we choose to report only 



3 we note that the secondary stars in CVs might be slightly ex- 
panded compared to those in detached systems. 
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Table 3. MORO Light Curve Fitting Parameters 



Fixed Parameters 



Parameter 


Solution #1 


Solution #2 


Solution #3 


Description 




(SOAK-BJ 


(Prompt-B) 


(Frompt-B VK1J 




Ti [K] 


28900 


28900 




28900 


primary effective temperature from spectroscopy 


Ai 


1.0 


1.0 




1.0 


primary bolometric albedo 


82 


0.0 


0.0 




0.0 


secondary radiation pressure parameter 


XlB 


0.26 


0.26 




0.26 


primary linear limb darkening coefficient, B filter 


XIV 








0.220 


primary linear limb darkening coefficient, V filter 


XlR 








0.190 


primary linear limb darkening coefficient, R filter 


x u 








0.165 


primary linear limb darkening coefficient, I filter 


91 


1.0 


1.0 




1.0 


primary gravity darkening exponent 


92 


0.32 


0.32 




0.32 


secondary gravity darkening exponent 


h 


0.0 


0.0 




0.0 


third light contribution 


A B [nm] 


435 


435 




435 


isophotal filter wavelength, B filter 


Av [nm] 








555 


isophotal filter wavelength, V filter 


A R [nm] 








640 


isophotal filter wavelength, R filter 


Ai [nm] 








790 


isophotal filter wavelength, I filter 


Free Parameters 


Parameter 


Solution #1 


Solution 


#2" 


Solution #3 a 


Description 




(SOAR-B) 


(Prompt- 


■B) 


(Prompt-BVRI) 




Q 


0.25 ± 0.04 


0.253 




0.250 


mass ratio 




79.75 ± 0.13 


79.61 




79.27 


inclination angle 


T 2 [K] 


2900 ± 500 


2861 




2724 


secondary effective temperature 


A 2 


1.90 ± 0.18 


1.84 




1.72 


secondary bolometric albedo 


Hi 


5.05 ± 0.09 


4.9937 




4.9701 


primary modified equipotcntial 


ft2 


2.85 ± 0.19 


2.8457 




2.7879 


secondary modified equipotcntial 


Si 


0.017 ± 0.010 


0.01971 




0.02152 


primary radiation pressure parameter 


Lib 


0.9999 ± 0.0011 


0.99998 




0.999992 


primary luminosity fraction 6 , B filter 










0.99993 


primary luminosity fraction, V filter 


Lin 








0.99979 


primary luminosity fraction, R filter 


L11 








0.99923 


primary luminosity fraction, I filter 


L2B 


0.00002 ± 0.00033 


0.000016 




0.000008 


secondary luminosity fraction 2 , B filter 


L2V 








0.00007 


secondary luminosity fraction, V filter 


L2R 








0.00021 


secondary luminosity fraction, R filter 


L21 








0.00077 


secondary luminosity fraction, I filter 


%2B 


0.601 ± 0.070 


0.734 




0.655 


secondary limb-darkening coefficient, B filter 


X2V 








0.866 


secondary limb-darkening coefficient, V filter 


%2R 








0.998 


secondary limb-darkening coefficient, R filter 










0.997 


secondary limb-darkening coefficient, I filter 


nja 


0.2058 ± 0.0013 


0.20782 




0.20833 


primary radius as fraction of orbital separation 


r-2/a 


0.1739 ± 0.0018 


0.17407 




0.17924 


secondary radius as fraction of orbital separation 




0.0028 


0.0121 




0.0160 


standard deviation about model fit 



a errors are not given for Solutions 2 & 3 since results are only derived using Solution 1. 
6 fraction of total system luminosity emitted by primary, defined as L — ' 
fraction of total system luminosity emitted by secondary, defined as — . 



conservative limits on the primary (sdB) and secondary (M 
dwarf) masses: Mi=0. 35-0.85 M and M 2 =0.11-0.17 M . 

Although we cannot claim any particular solution as the 
definitive one, we continue our analysis and discussion using 
the best-fitting MORO solution (the red point in Figure 4) 
to illustrate a set of self-consistent parameters that might 
describe EC 10246. For completeness, Table 3 presents the 
fixed and adjustable parameters corresponding to the best- 
fitting MORO solutions for all data sets. The system param- 
eters of the best-fitting models to each of the three data sets 
are compatible with one another. Throughout the rest of the 
manuscript, we use only Solution #1 (the highest-quality 



solution) for illustration purposes. Error bars are shown in 
Table 3 for this solution only and represent the standard 
deviations of the best-fitting parameters from the family of 
27 models MORO found for the SOAR light curve. 

Figure 5 presents the best-fitting model solution to the 
SOAR light curve, which has a mass ratio q — 0.254 and 
nicely fits the data. We also show the result of the simul- 
taneous PROMPT BVRI fit in Figure 6. Of all the output 
parameters, the relative radii (n, r 2 ) and inclination an- 
gle (i = 79°. 75) are the most tightly-constrained given their 
strong dependence on the relatively-easy-to-measure eclipse 
depths and durations. Figure 7 shows a scaled model of the 
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best-fit solution over one orbital cycle in 0.125-phase incre- 
ments, starting at the secondary eclipse. The large reflec- 
tion effect can only be modeled using a secondary albedo 
A2 — 1.9, a physically unrealistic value. Other attempts to 
model sdB+dM light curves have run into similar problems, 
which probably illustrates a mishandling of the re-radiated 
light in the theoretical models; to compensate for this short- 
coming and still fit the light curve accurately, the secondary 
albedo can be adjusted. For et al. (2010), for instance, find 
acceptable fits to 2M 1533+3759 only when the secondary 
albedo is adjusted to A2 = 2.0. The frequency redistribution 
in the M dwarf photosphere, which leads to emission primar- 
ily in the B filter bandpass, might produce this effect. We 
discuss this further in §5. 

Finally, after subtracting Solution # 1 from the SOAR 
light curve, we computed the Fourier transform of the resid- 
uals to look for rapid pulsations in the sdB star. We find 
no significant signals in the Lomb-Scargle periodogram and 
limit the amplitudes of rapid, coherent photometric oscilla- 
tions to < 0.08%. 



4 SYSTEM PARAMETERS 

The mass function for EC 10246, defined by, 
K'(P Mig 3 sin 3 i 



f 



2tvG 



(i + g) 2 



(2) 



has a value of 0.0045 ± 0.0003 Mq, as calculated from the 
period and sdB velocity. After inserting into the above equa- 
tion the mass ratio and inclination angle from the best- 
fitting MORO model, we find Mi= 0.45 ± 0.17 Mq and 
Mi = 0.12 ± 0.05 Mq. The orbital separation, then, must 
be a — 0.84 ± 0.10 Rq, according to Kepler's Third Law. 
With this value, we can convert the radii from the MORO 
output, which were given in fractional units of a, to phys- 
ical units. We find i?i = 0.17 ± 0.02 R Q and R 2 = 0.146 
± 0.018 Rq, consistent with an sdB and M dwarf. Under 
the assumption of synchronous rotation, we calculate a ro- 
tational velocity of Kot = 74 ± 9 km s _1 for the sdB com- 
ponent using the above orbital parameters. Such a rotation 
would give rise to a 1- A broadening of the subdwarf spectral 
features that is four times smaller than the spectral resolu- 
tion of our SOAR/Goodman spectra and, thus, difficult to 
measure using the current data set. Table 4 summarizes the 
derived system properties for EC 10246. 

We stress that T2 is poorly constrained since the sec- 
ondary contributes a negligible amount of the system light 
in the optical (aside from the reflection effect). Nonetheless, 
T2, M2, and R2 are roughly consistent with the theoreti- 
cal models of single lower-main sequence stars derived by 
Baraffe et al. (1998), in addition to the empirical spectral 
type-radius relation determined by Rebassa-Mansergas et 
al. (2007) for M-dwarfs in post common envelope binaries. 
Owing to the relatively large uncertainties in M2, R2, and 
T2, it is difficult to determine the precise spectral type of the 
cool companion in EC 10246-2707; both theoretical and em- 
pirical M-R relations point to a companion near or slightly 
later than M5-6V. The secondary's mass and radius are also 
consistent with the empirical M-R relation for CV secon- 
daries from Knigge et al. (2011); we note, however, that the 



Table 4. System Properties from the Best— Fitting Solution 



System Properties 



p 


0.118 507 993 6 


± 0.000 000 000 9 


days 


To 


2 455 680.562160 


± 0.000 016 


BJD 


q 


0.25 


± 0.05 




a 


0.84 


± 0.10 


R 


i 


79.75 


± 0.13 


deg 


d 


830° 


± 50 


pc 



Primary Properties 

(subdwarf B star) 



Mi 

T e ff 
log g 

lQ S N(H) 

Ki 

Vrot 



0.45 

0.17 

28900 6 

5.64 6 

-2.46 6 

71.6 6 

74 c 



± 0.17 
± 0.02 
± 400 
± 0.02 
± 0.15 
± 1.7 
± 9 



R 
K 



km s 
km s 



Secondary Properties 

(M4-M5 dwarf) 



M 2 


0.12 


± 0.05 


M 


R2 


0.146 


± 0.018 


R© 


T e g 


2900 


± 500 


K 


log g 


5.2 d 


± 0.3 


cm s~ 2 


K 2 


286 e 


± 5 


km s _1 


Vrot 


63 c 


± 8 


km s _1 


A 2 


1.90 


± 0.18 





a estimated from T\ , Ri , and extinction (Schlcgel et al. 1998) 
b measured from spectroscopy (not from light curve modeling) 
c assumes rotation synchronised with orbit 
d calculated from M2, R2 
e calculated from q and K\ 



M dwarfs in such systems might be slightly expanded com- 
pared to detached systems such as EC 10246-2707. Lastly, 
we draw attention to the close coincidence of our secondary 
mass and radius with those of Proxima Centauri (GJ 551) 
measured by Demory et al. (2009) using the VLTI interfer- 
ometer and VINCI; they classify the star as an M5.5 dwarf. 



5 DISCUSSION 

In Table 5, we list the parameters of all known eclipsing 
sdB+dM systems, ordered by increasing orbital period. Our 
preliminary light curve modeling solution implies EC 10246 
is a typical sdB+dM system. Two-thirds of the M dwarfs 
in eclipsing sdB+dM systems have reported spectral types 
of M4-M5, according to their light-curve modeling solu- 
tions. Only NSVS 14256825 (M2.5-M3.5; not well-studied), 
J08205+0008 (>M10), and HS 2231+2441 (brown dwarf; 
not well-studied) differ in this respect. Whether all eclips- 
ing sdB+dM binaries contain ~ M4-M5 dwarfs remains to 
be seen. It is interesting to note, however, that these spec- 
tral types correspond to the boundary at which M dwarfs 
become fully-convective. 

Although our study is not the first to find a sec- 
ondary albedo A2 > 1.0 (For et al. 2010), the value we 
find from the best-fitting MORO solution is still discon- 
certing. We can leverage HW Vir's similarity to EC 10246 
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Figure 7. Snapshots from the best-fitting MORO model for EC 10246-2707 over one orbital cycle in 0.125-phase increments starting 
from the secondary eclipse (top left), progressing through the primary eclipse (centre image), and returning to the secondary eclipse 
(bottom-right; sequence progresses from left-to-right, top-to-bottom). The system is viewed from an inclination angle of 79.8°. The 
separation distance and stellar radii are drawn to scale. 



ct> = 0.25 



$ = 0.375 



( 



to perform a consistency check by comparing their reflec- 
tion effect amplitudes, similar to what For et al. (2010) did 
for 2M 1533+3759. We estimate the reflection effect ampli- 
tude using a simplified model in which the orbital plane is 
viewed exactly edge-on, both stars are perfect spheres, the 
secondary emits no light of its own, and the primary uni- 
formly illuminates exactly half of the secondary's surface. 
In this model, the reflection amplitude (given as a fraction 
of the primary flux over some wavelength range) scales with 
A2R2/C1 2 . Under the assumption that the reflection effect 
physics is identical for two binary systems (j, k), the ratio 
of their secondary albedos will be given by the expression 




From our estimates of R2, and a, and Lee et al.'s (2009) 
values for HW Vir, we expect a reflection amplitude 70% 
that of HW Vir in the B-filter. The observed amplitude is 
93% of its value, implying the reflection effect actually is 
stronger in EC 10246 and requires a larger albedo to model 
properly. Of course, the value of R2 used to demonstrate this 
inconsistency is somewhat dependent upon A2, so this line 
of reasoning should be taken likely. A scaling comparison to 
2M 1938+4603 (0stensen et al. 2010) gives a similar result. 

Eclipsing sdB+dM binaries are expected to merge even- 
tually due to orbital angular momentum loss, primarily from 
the emission of gravitational radiation and magnetic brak- 
ing from the M dwarf. Once their orbits shrink sufficiently 
so that the secondary fills its Roche lobe, mass transfer will 
begin. If the cool companion in EC 10246 truly is of spec- 
tral type M5V or later, then it is fully convective and mag- 



netic braking should be relatively inefficient (Schreiber et al. 
2010). Consequently, gravitational radiation will dominate 
the rate of angular momentum loss, and the system should 
initiate mass transfer long after the sdB has evolved into a 
white dwarf and become a cataclysmic variable (Schreiber 
& Gansicke 2003). We can estimate the orbital period at 
this transition point (-Pcv) using Kepler's law. The ratio of 
the current period to P cv equals the ratio of R\f to -R2' 5 , 
where J?2L is the current Roche lobe radius of the secondary 
(Eggleton 1983). We calculate a period of P cv ^ 1.48 firs, 
which falls above the minimum for CVs (~ 80 min) but 
below the lower-edge of the famous period gap (~ 2 hrs). 

Precise eclipse measurements over a long timespan can 
constrain the orbital decay rates of sdB+dM binaries. Esti- 
mates of P have been reported for many systems, including 
NSVS 14256825 (+1.09 x 10 -10 s s _ \ Kilkenny & Koen 
2012 4 ), HW Vir (-2.3 x f0~ n s s" 1 , Lee et al. 2009), and 
NY Vir (-IT x f0" 10 s s _1 , Kilkenny 201f ; Qamurdan et al. 
2012). For all of these systems, however, it is probably too 
early to know whether the parabolic terms in their O-C di- 
agrams are instead long-term cyclic variations. Even with 
15 years of eclipse timings for EC 10246, we are not able 
to claim a statistically significant P at this time. Given the 
substantial amount of scatter in our O-C diagram, addi- 
tional high-precision timings (e.g., Parsons et al. 2010) are 
needed to rule out small-amplitude changes, such as those 
claimed for NSVS 14256825 (Kilkenny & Koen 2012; Beuer- 
mann et al. 2012). We note that AA Dor also shows no de- 



4 Beuermann et al. (2012) claim a third body in a highly elliptical 
orbit is responsible for the observed O-C changes 
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tectable period change after more than 35 years of monitor- 
ing (Kilkenny 2011). If the physics governing orbital shrink- 
age is the same in all HW Vir-type systems, then the ap- 
parent orbital stability of AA Dor and EC 10246 suggests 
that the parabolic trends claimed in sdB+dM O-C diagrams 
might instead be combinations of one or more long-term 
cyclic changes. 

Finally, we draw attention to the opportunity eclipsing 
sdB+dM systems provide to measure stellar masses in a rel- 
atively model-independent way. Due to the finite speed of 
light, it is possible to calculate both masses in an eclipsing 
binary (if Mi 7^ M2) by measuring the time delay of the 
secondary eclipse with respect to one half-period after the 
primary, a delay resulting from the extra light travel time 
(see Kaplan 2010). For a circular orbit, the delay is given by 

A* LTT = ^i(I-l) (4) 

and increases with decreasing mass ratio. This technique 
works for eccentric systems, too, but since non-zero eccen- 
tricity also affects the relative spacing between eclipses, the 
orbital geometry must be determined first with high pre- 
cision. The mass ratio can be computed independently of 
the inclination angle, while adding an estimate of this value 
(which is well constrained in eclipsing systems) allows one 
to compute the individual masses. Barlow et al. (2012) and 
Bloemen et al. (2012) recently applied this method to 2M 
1938+4603 (sdB+dM) and KOI- 74 (A dwarf+ white dwarf), 
respectively, using Kepler eclipse timings to derive the com- 
ponent masses. With their relatively small mass ratios and 
assumed circular orbits, eclipsing sdB+dM systems are good 
candidates for this technique. The last column of Table 5 
shows the predicted time delays for the known eclipsing 
sdB+dM systems. Most fall near 2 s and are observable given 
an adequate number of precise eclipse timings (as long as the 
eccentricity is well-constrained). Measurements of masses 
derived in this way can help test the outputs of binary light 
curve-modeling solutions and, in systems with pulsating sdB 
components (like PG 1338-018 and 2M 1938+4603), astero- 
seismology. 



6 SUMMARY 

We have presented photometric and spectroscopic observa- 
tions of EC 10246-2707, a previously-unpublished eclipsing 
sdB+dM binary. We find an orbital period of 2.84 hr and 
an sdB orbital velocity of 71.6 km s" 1 , which together give 
a mass function equal to 0.0045 Mq. Parameters derived 
from modeling the light-curve with MORO imply the hot 
subdwarf has a mass near the canonical value, while the 
cool companion is likely an M dwarf with spectral type near 
M5V. The best-fitting model gives sdB and M dwarf masses 
of 0.45 Mq and 0.12 Mq, respectively. Due to numerous 
assumptions made during the modeling and the degener- 
acy associated with these light-curve solutions, however, we 
cannot claim a unique solution at this time. 

Like other HW Vir-type systems, the orbit of EC 10246 
will decay over time, and eventually the system will be- 
come a cataclysmic variable with an orbital period below 
the CV period gap. Long-term monitoring of the eclipse tim- 



ings might, in due course, provide constraints on the rate of 
period change, as long as the effects of the system's proper 
motion can be identified and removed if necessary. Timing 
measurements also have the potential to uncover the pres- 
ence of tertiary members in the system down to planetary- 
size, as reported for PG 1336-018 (Qian et al. 2011), HW Vir 
(Lee et al. 2009), NSVS 14256825 (Beuermann et al. 2012) 
and HS 0705+6700 (Beuermann et al. 2012; Qamurdan et al. 
2012). 
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Table 5. System Parameters of Known sdB-f-dM Eclipsing Binaries, with sdB Properties 



System Name 


P 


Ki 


a 


i 


Mi 


M 2 




log g 


A*ltt 


Reference 


(sorted by P) 


[hr] 


[km s" 1 ] 


[R©] 


[deg] 


[M ] 


[M ] 


[K] 


[cm s~ 2 ] 


[s] 




J162256+473051 a 


1.8 


- 


- 


- 


- 


- 


- 


- 


- 


Schaffenroth et al. (2011) 


HS 0705+6700 


2.29 


85.8 


0.81 


84.4 


0.48 


0.13 


28800 


5.40 


2.0 


Drechsel et al. (2001) 


J08205+0008 b 


2.30 


47.4 


0.6/1.1 


85.90 


0.25/0.47 


0.045/0.068 


26700 


5.48 


2.5 


Geier et al. (2011) 


PG 1336-018 


2.42 


78.6 


0.723 


80.67 


0.389 


0.110 


31300 


5.74 


1.8 


Vuckovic et al. (2007) 


NSVS 14256825 


2.65 


73.4 


0.74 


82.5 


0.35 


0.097 


42300 


5.50 


3.4 


Almeida et al. (2012) 


HS 2231+2441 


2.65 


49.1 


1.18 


79.1 


0.265 d 


0.05 d 


28370 


5.39 


2.9 


0stensen et al. (2007) 


HW Vir 


2.80 


82.3 


0.8594 


80.98 


0.485 


0.142 


28488 


5.60 


2.1 


Lee et al. (2009) 


EC 10246-2707° 


2.86 


71.6 


0.84 


~80 


~0.45 


~0.12 


28900 


5.64 


2.3 


this manuscript 


BUL-SC16 335 


3.00 


















Polubek et al. (2007) 


2M 1938+4603 


3.02 


65.7 


0.89 


69.45 


0.48 


0.12 


29564 


5.43 


2.3 


0stensen et al. (2010) 


ASAS 10232-3737.0 C 


3.34 


81.0 


0.963 


65.86 


0.461 


0.157 


28400 


5.56 


4.1 


Schaffenroth et al. (2011) 


2M 1533+3759 


3.88 


71.1 


0.98 


86.6 


0.376 


0.113 


29230 


5.58 


2.5 


For et al. (2010) 


AA Dor" 


6.26 


39.2 


1.28 


89.21 


0.47 


0.079 


42000 


5.46 


7.2 


Klcpp & Rauch (2011) 



a unclear whether M dwarf or brown dwarf companion. 
b likely a brown dwarf, not an M dwarf. 
c marginal eclipse, large solution errors. 

d For et al. (2010) argue that the log g cited by 0stensen et al. (2007) is most likely underestimated, in which case both the sdB and 
companion are more massive than reported. 

e best-fitting solution reported; actual parameters may vary significantly. 
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